On Camassa-Holm equation with self-consistent sources and its 

solutions 



Yehui Huang* Yuqin Yao^and Yunbo Zeng^ 
Department of Mathematical Sciences, 
Tsinghua University, Beijing, 100084, P.R. China 



Abstract 

Regarded as the integrable generalization of Camassa-Holm (CH) equation, the CH equa- 
tion with self-consistent sources (CHESCS) is derived. The Lax representation of the CHESCS is 
presented. The conservation laws for CHESCS are constructed. The peakon solution, N-soliton, 
N-cuspon, N-positon and N-negaton solutions of CHESCS are obtained by using Darboux trans- 
formation and the method of variation of constants. 
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1 Introduction 

Camassa-Holm (CH) equation, which was implicitly contained in the class of multi-Hamiltonian 
system introduced by Fuchssteiner and Fokas ^ and explicitly derived as a shallow water wave 
equation by Camassa and Holm has the form 

Ut -\- 2uJUx — Uxxt + 3uUx = lUxUxx UUxxxi (1.1) 

where u = u{x,t) is the fluid velocity in the x direction and the constant 2uj is related to the critical 
shallow water wave speed. Let q = u — Uxx + w, we have the following equivalent equation ^. 

qt + 2uxq + uqx = 0. (1.2) 

It was shown by Camassa and Holm that this equation shares most of the properties of the integrable 
system of KdV type It possesses Lax pair formalism and the bi-hamiltonian structure. When 
w > 0, the CH equation has smooth solitary wave solutions. When w — > 0, these solutions become 
piecewise smooth and have cusps at their peaks. These kind of solutions are weak solutions of (jl.2p 
with uj = and are called "peakons". Since the works of Camassa and Holm, this equation has 
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become a well-known example of integrable systems and has been studied from many kinds of views 

4-12 

Soliton equations with self-consistent sources (SESCS) have attracted much attention in recent 
years. They are important integrable models in many fields of physics, such as hydrodynamics, 
state physics, plasma physics, etc ^3-25^ ^or example, the KdV equation with self-consistent sources 
describes the interaction of long and short capillary-gravity waves The nonlinear Schrodinger 
equation with self-consistent sources represents the nonlinear interaction of an electrostatic high- 
frequency wave with the ion acoustic wave in a two component homogeneous plasma The KP 
equation with self-consistent sources describes the interaction of a long wave with a short wave 
packet propagating on the x-y plane at some angle to each other The SESCS were firstly 
studied by Melnikov A systematic way to construct the soliton equations with self-consistent 

sources and their zero-curvature representations is proposed ^1-24 fpj^g problem of finding soliton 
solutions or other specific solutions for SESCS has been considered in the past by many authors 

13-25 

The present paper falls in that line of the work on the CH equation concerning with establish- 
ing the many facts of its completely integrable character, aiming at the integrable generalization 
of CH equation by deriving the Camassa-Holm equation with self-consistent sources (CHESCS) 
and finding its solutions. We first construct the CHESCS by using the approach presented in the 
reference ^1-24^ The Lax pair of the CHESCS is obtained, which means that the CHESCS is Lax 
integrable and can be viewed as integrable generalization of CH equation. Since the CH equation 
describes shallow water wave and the SESCSs in general describe the interaction of different soli- 
tary waves, it is reasonable to speculate on the potential application of CHESCS, that is, CHESCS 
may describe the interaction of different solitary waves in shallow water. It was pointed out 
that SESCS can be regarded as soliton equations with non-homogeneous terms, and accordingly 
proposed to look for explicit solutions by using the method of variation of constants. Applying 
this technique to CHESCS we have been able to find its peakon solution. In order to find other 
solutions of CHESCS, we consider the reciprocal transformation ^8:29^ which relates CH equation to 
an alternative of the associated Camassa-Holm (ACH) equation, and propose the reciprocal trans- 
formation, which relates the CHESCS to associated CHESCS (ACHESCS). By using the Darboux 
transformation (DT), one can find the n-soliton and n-cuspon solution ^'^ as well as positon and 
negaton solution of alternative ACH equation. Then by means of the method of variation of con- 
stants, we can obtain the N-soliton, N-cuspon, N-positon and N-negaton solution for ACHESCS. 
Finally, using the inverse reciprocal transformation, we obtain the N-soliton, N-cuspon, N-positon 
and N-negaton solution of CHESCS. 

This paper is organized as follows. In section 2, we present how to derive the CHESCS and 
its Lax representation. In section 3, the conservation laws of the CHESCS are constructed. In 
section 4, the peakon solution is obtained. In section 5, we consider the reciprocal transformation 
for CH equation and CHESCS, respectively. In section 6, by using the DT, we find the solution for 
alternative ACH equation, then by using the method of variation of constants and inverse reciprocal 
transformation, we obtain the N-soliton, N-cuspon, N-positon and N-negaton solution for CHESCS. 
In section 7, the conclusion is presented. 



2 



2 The CHESCS and its Lax pair 



2.1 The CHESCS 

The Lax pair for CH equation (II. 2p is given by ^ 



^2A 2 

It is not difficult to find that 



^xx = (Ag + ^)v', (2.1a) 
1 1 

'^t = (ttt - '^)'^x + -Ux^- (2.1b) 



6q 

The CH equation possesses bi-hamiltonian structure ^ 



^ = -A(^^ (2.2) 



6 Ho ^5 Hi 



where 



Ho = / ti^ + n^dx, 



Hi = - I u'^ + uu'tdx. 



K = -d^ + d, 

J = dq + qd, 
1 
2 
1 
2 

According to the approach proposed in the reference ^1-24^ CHESCS is defined as follows 

N 

= -{qd + dq){u + 2Y,>^jV]) 
i=i 

N 

= -2qUx-uqx + '^{-8Xjqipjipjx-2Xjqxipj), (2.4a) 

^j,xx = {Xjq + j = 1, • • • ,N, (2.4b) 
which has a equivalent form by using (2.4b) 

N 

qt = -2qux - uqx + '^[i(p'j)x - {^'j)xxx], (2.5a) 

^j,xx = {Xjq + -)(pj, j = l,---,N, (2.5b) 
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2.2 The Lax representation of the CHESCS 

Based on the Lax pair of the CH equation (2.1), we may assume the Lax representation of the 
CHESCS (2.4) or (2.5) has the form 

^xx = {M + (2-6a) 
ift = -^Br,(p + (2.6b) 

j=l 3 j=l 

where f{(pj) is undetermined function of ipj. The compatibihty condition of (2.6a) and (2.6b) gives 

Xqt = LB + X{2B^q + Bq^), (2.7) 
where L = -^d^ + ^d. Then (2.6) and I^Tii yields 

^<it = -^ E T^[/"Vix + 3(/'V, - /')(A,9 + j)^,x + XMf'^j - 2/)] 



A - A,- 

j=l 



N N 



1 ^ 

^ + + Xjflx^j - f'^j] + '^ajiQxf + 2qf'ipj^). (2. 



Here /' denotes the partial derivative of the function / with respect to the variable ipj. In order 
to determine /, aj and (3j, we compare the coefficients of j^x', A and A", respectively. We first 

observe the coefficients of , then the coefficients of ip^^, (pjx and other terms gives rise to, 
respectively 

/"'=0, /V^.-/' = 0, /V,-2/ = 0, 
which leads to / = Substituting / = bip'^ into the coefficients of A in (2.8) gives 

N N 
qt = -2qux - uqx + 4g E Pj^^j^jx + Qx E 1^3^^]- 

3=1 j=l 

Comparing the above equation and (2.4a), we can determine 

b = -2, Pj = Xj. 

Substituting / = —2f'j, and Pj = Xj into the coefficients of A'' in (2.8), we obtain 

aj = X]. 
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Thus we obtain the Lax pair of the CHESCS (2.5) 



'Pxx = + M)'~P^ (2-9a) 

1 ^ AA- • 

'ft = Y</0 + ^2A ~ ^^"^"^ + ^Xl X^A^^'^-''^'^ ~ 'fd'P^^)- (2-9b) 

i=i 



which means that the CHESCS (2.5) is Lax integrable. 



3 The infinite conservation laws of the CHESCS 

With the help of the Lax representation of the CHESCS, we could find the conservation laws for 
the CHESCS by a well-known method. First we assume that u, and its derivatives tend to 
when |x| — > oo. Set 

r = — , (3.1) 

then the identity 

d dXvup d din if 

together with (2.10) implies that CHESCS has the following conservation law: 

d d .iff. d ,1 XXj ,,1 , x"-^ AA,- 9,^, , , 

. Using (2.10a) gives rise to 

T:, = ^ + qX-T^. (3.3) 

Let 

oo 

v=Y,i^m>r^, (3.4) 

m=0 

then /Xto is the density of conservation laws. 
Define 

1 ^ AA 1 ^ AA °° 1- 

2«x + 2^ ^^^j^j. + ((^ - u) - 2^ ^^<^2)p ^ ^ ^^^i-H^ (3 5) 

j=l j=l m=0 

It is foTind that the density of the conservation laws jim aiid the flux of the conservation laws 
Frfn satisfy the following recursion relation: 

M = Vq, 

1 Qx 

m = -7 — . 

4 q 

^-vm— 1 

2/xo 
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N 

N I ^ 

m TV 
j=0 j=l 

m N N 

F2m+1 = ^(^i»+2j^A}+V')/"2m-2m+2 5^AJ^+Vi¥'jx, ^ > 1, (3.7) 
j=0 i=l j=l 

where u^'^^ = u, u^^^ = 1, u^*) =0, i > 1. 

After some calculations we can find the first few conserved quantities given by /xq, /X2 and ^4 
are as follows 



H_i = j y/^dx, (3.8a) 
The corresponding flux of the conservation laws are 

N 

G_i = {-u-2Y,^jv])VQ, (3-9a) 

i=i 

G_2 = (l + 2E\V,^)V^+(^ + 2EV,')(l^(7= + -|^)-(^U (3.9b) 

+ I^(1 + 2EaM)(;7^ + ^) + 2EaM^/^. (3.9c) 

As the space part of the Lax Pair of the CHESCS is the same as that of CH equation, the 
densities of the conservation laws of the CHESCS are the same as those of the Camassa-Holm 
equation As the time part of the Lax pair is different, the fluxs of the conservation laws for CH 
equation and CHESCS are different. 

4 One peakon solution of the CHESCS 

The CH equation (|1.2p has peakon solutions ^ 

u = ce-l^-'=*+°l, (4.1) 
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where a is an arbitrary constant. The corresponding eigenfunction of (2.1) is 

(^ = /Je-5l^''^*+"l, (4.2) 

where (3 is an arbitrary constant. 

Since the CHESCS (2.5) can be considered as the CH equation ()1.2p with non-homogeneous 
terms, we may use the method of variation of constants to find the peakon solution of CHESCS 
from the peakon solution (4.1) and (4.2). Taking a and /? in (j4.ip and (4.2) to be time-dependent 
a{t) and /3(t) and requiring that 

u = ce-'^-'^*+"WI, (4.3a) 
ip = /3(t)e-5l^-^*+"WI (4.3b) 

satisfy the CHESCS (2.5) for iV = 1. We find that c = i, a{t) can be an arbitrary function of t 
and P{t) = yja'{t)c. So we have the one peakon solution for (2.4) with = 1, Ai = A = ^ 

^ = v/a^e-5l^-^*+°(*)l (4.4b) 

The one peakon of the CHESCS also has a cusp at its peak, located at x = ct — a{t). We note 
that for the one peakon solution of the CH equation, the solution travels with speed c and has a 
cusp at its peak of height c, for the CHESCS, the cusp is still at its peak of height c, but the speed 
c — of the wave is no longer a constant. 

5 A reciprocal transformation for the CHESCS 

Let r = yfq, by the reciprocal transformation ^>28.29 

dy = rdx — urds, ds = dt, 
and denoting / = r^^cj), the Lax pair (2.1) of CH equation is transformed to the following system 

(^yy = {X + Q + ^)^, (5.1a) 

where 

Q = -i(!l)2 + ^ + JL_ J_. (5.2) 
4 r 2r 4r2 4u; 

The compatibility condition of (5.1a) and (5.1b) gives an alternative of the associated CH (ACH) 
equation 

Qs — fy-i (5.3a) 
111 

- -^ry + -ryyy " -QyV " QVy = 0. (5.3b) 
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We now consider the reciprocal transformation for the CHESCS (2.5). (2.5a) gives 

N 

rt = -{ru),-2Y,^j{ry^%- (5-4) 

(j5.4|l shows that the 1-form 

N 

oj = rdx — {ru + 2 Xjrip'j)dt (5-5) 
i=i 

is cfosed, so we can define a reciprocal transformation (x, t) (y, s) by the relation 

N 

dy = rdx — {ru + 2 Xjrip'j)ds, ds = dt, (5-6) 

i=i 

and we have 

Denoting ip = r~'2^^ ipj = r~2ipj and using (j5.2p . the Lax pair (2.9) of CHESCS (2.5) is corre- 
spondingly rewritten as 



V'yj/ = (A + g + ^)'0, (5.8a) 



N -.2 



1 1 v-^ ^jVi 

= l^ir^y - ^^s/V') + 2 Y3T-(V'j?;V' - ^jV's/)- (5.8b) 



2A' 2 ^ A- A,- 

The compatibility condition of (5.8a) and (5.8b) leads to an associated CHESCS (ACHESCS) 

N 

Qs = ry - X'jipjipjy, (5.9a) 
i=i 

111 

- -^ry + -ryyy " -Q^F " Qr^, = 0, (5.9b) 

^jyy = (Aj +Q + i = 1, 2, . . . , iV. (5.9c) 



The Eqs. (j5.9p can be regarded as the Eqs. (|5.3p with self-consistent sources. In order to obtain the 
solutions of the CHESCS (2.5), we have to get the relation of the variables (y, s) and the variables 
{x,t). From the reciprocal transformation, we have 

By making use of the compatibility of the above two equations, we have 

xiy,s) = I ^dy. (5.11) 
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The solutions of the CHESCS (2.5) with respect to the variables (y,s) are given by 

■tp ■ 

q = r'^{y,s), (fj{y,s) = (5.12a) 

TV N 

u{y, s)=r^-rys + ^- 2rr, ^ Xj{^])y - 2r^ Yl ^M)yy " ^' (5-12b) 

xiy,s)=J^dy. (5.12c) 

We now prove (5.12b). From q = u — Uxx + ^ and the reciprocal transformation ()5.7p . we have 

u = q + rryUy + r'^Uyy — uj. (5.13) 
By using the reciprocal transformation ()5.7p . ()5.4p gives rise to 

N 

i=i 

Substituting (|5l4ll and ([52]) into (f5l3]) leads to (5.12b). 

6 The solutions for the CHESCS 

Notice that Q = 0, r = ^/uj is the solution of (5.2) and (5.3). Let the functions (j)Q{y,s,X), 
^'i(y, s,Ai), ••• ,^'„(y,s,A„) be different solutions of (5.1) with Q = 0, r = and the corre- 
sponding A and A = Ai, • • • , A„, respectively. We construct two Wronskian determinants from these 
functions 

W, = W{^^,--- Mr\^2r-- ,*n,--- MrT"^), (6.1a) 

W2 = W{^U--- ,^^r\'^2,--- '^n,--- ,^^r"\'/'0), (6.1b) 

where mj > are given numbers and ^'j*'' = ^-^^jT^^\x=Xy 

Based on the generalized Darboux transformation for KdV hierarchy and using (5.3a), the 
following generalized Darboux transformation of (5.1) is valid ^'30,31 



Q{y,s) = -2^logWi, (6.2a) 



r{y,s) = (6.2b) 
4>{y,s,X) = -^, (6.2c) 



namely Q{y,s), r{y,s) and (p{y,s,X) satisfy (5.1), (5.2) and (5.3) 
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6.1 The multisoliton solutions 



Take and $j be the solutions of Eq.(5.1) with Q = 0, r = ^/uj and \ = — < 0, or 
Aujkf - 1 < 0, (0 < fci < A;2 < • • • < ^n) as follows 

= cosh^i, i is an odd number, (6.3a) 

= sinh^i, i is an even number. (6.3b) 

$i = (6.4) 

where henceforth 

Ci = h[y+ ^^^2 _ I + • (6-^) 

By using Darboux transformation (6.2) with mi = • • • = m„ = 0, the n-soliton solution Q{y, s) and 

r{y,s) of (5.3) and the corresponding eigenfunction 4>i{y,s,Xi) of (5.1) with Aj = kf — is given 

by 

Q{y, s) = -2[logW{^i,^2, • • • , ^n)]yy, (6.6a) 

r{y, s) = - 2[logW{^i, ^2, • • • , ^n)]ys, (6.6b) 



When n = 1 and 4kfuj-l < 0, (6.6) gives rise to one soliton solution for (5.3) and the corresponding 

4aj 



eigenfunction of (5.1) with \i = kf — 



Qiy, s) = -2klsech^ii, r{y, s) = ^/uj - ^^1^ ^^^^ ^1 ^ 7^^) 

Akfuj — 1 

= kisech^^i. (6.7b) 

Since Eq. (j5.9p can be considered to be Eq. (|5.3p with non- homogeneous terms and (pi satisfies 
(5.1a) with A = Ai, we may apply the method of variation of constant to find the solutions of the 
CHESCS (j5.9p by using the solution ()6.7p of ACH equation (j5.3p and corresponding eigenfunction. 
Taking ai in (6.5) to be time-dependent functions ai{s) and requiring that 

Q{y,s) = -Iklsech^ii, f{y,s) = y/LO ^-^ ; — , (6.8a) 

Ak{uj — 1 

V'l = I3i{s)kisech£^i (6.8b) 
satisfy the system ()5.9p for = 1, henceforth, we denote 

C^ = h[y+ -^2 — 7 + «'(^)]- (6-9) 

^UJK^ — 1 

We find that ai{s) can be an arbitrary function of s and 



Pi{s) = Y^^V2«'i(^)- (6.10) 
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So the one-soliton solution of the CHESCS (2.5) with = 1 and Xi = kf — ^ < is obtained 
with respect to the variables {y, s) from (j5.12p 

q{y,s)=uj{l —2 — ) , (6.11a 

a; — 1 

Sk'^uj'^ sech'^^i 

^ (1 - 4kfoj){l - ikfuo + AkfoJsech^Ci) ' ^^'^^^^ 

2Y^2a'i(s)fcia;sec/i^i 1 1 \ 

9^i(?/,s)= / (6.11c) 
V^(l - 4A;2l^)(1 - 4A:2w + 4A:^wsec/i2^i) 



r \ y o^ 1 - 2A;iV'^tanh^i /ruja 
2;(y,s) = ^ -21n— — — - . (6.11d) 

V 1 + 2A;i y tanh 

The requirement Aklu — 1 < guarantees the nonsingularity of solution (6.11). 
In Fig 1, we plot the single soliton solution of u and ipi. 




Figure 1. Single soliton solutions for u and the eigenfunction ipi when w = 0.01, fci — 1, ai(s) — 
4s, s = 2. 

When n = 2, \i = kf - ^ < 0, X2 = k'^ - ^ < 0, we have 

^i = cosh^i, ^2 = sinh(,2, (6.12a) 
$1=6^1, ^>2 = e^^ (6.12b) 
H^i(^'i, ^'2) = k2Cosh(,2Cosh(,i — kisinh^2sinh^i, (6.12c) 
^2(^1, ^2, ^1) = k2{k^ - kj)sinh^u (6.12d) 
1^2(^1, ^2, ^2) = ki{kl - k1)cosh^2, (6.12e) 

(6.12f) 

Then (6.6) with n = 2 gives rise to two soliton solution for (5.3) and the corresponding eigenfunction 
of (5.1). In the same way as we did on the one-soliton solution, we can apply the method of variation 
of constants to get the two soliton solution of the ACHESCS (2.9) which together with (5.12) yields 
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to the two soliton solution for CHESCS (2.5) with N = 2, Xi = kj - ^, X2 = - ^ 

r{y,s) = V^- 2[Zo5Wi(*i,*2)]j/.|^,=^-,, (6.13a) 

= J^^RZ^Sa^l^l , = 1,2. (6.13b) 
(l-4fcfa;) /n(fc|-A;2)W^i(^'i,*2) 

In Fig 2 wc plot the interactions of two soliton solution for u and <^i, (^2, which is shown that u is 
elastic collision. 




Figure 2. two soliton solutions for u and the eigenfunction ipi, 1^2 when w = 0.01, fei =2, k2 = 1, q:i(s) = 2s, 0:2(5) 



Notice that the soliton solutions of CHESCS contains arbitrary s functions aj{s). This implies 
that the insertion of sources into the CH equation may cause the variation of the speed of soliton. 
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In the same way as in the reference we may apply the method of variation of constant to 
find the N-sohton solution of (2.5) with Xi = kf — > 0, i = 1, - ■ ■ ,N, from (|5.12p . where 

r{y,s) = VII] -2[logWi{^i,^2,--- , ^'jv)]ys|g,=^- , (6.14a) 
6.2 The multicuspon solutions 

Take and be the solutions of Eq. (|5.ip when Q = 0, r = y/uj and \ = ~ ^ > ^ 
(0 < ki < k2 < • ■ • < kn), as follows 

= sinh^i, i is an odd number, (6.15a) 
= cosh^i, i is an even number. (6.15b) 

$i = e«% (6.16) 

where is given by (6.5). 

The n-cuspon solution Q{y, s) and r(y, s) of (5.3) and the corresponding eigenfunction (j)i{y, s, Aj) 
of (5.1) with Aj = kf — is given by (6.6). 

When n = 1 and 4:kfuj — 1 > 0, (6.6) gives rise to one cuspon solution for (5.3) and the 
corresponding eigenfunction of (5.1) with = kf — 

Q{y, s) = 2klcsch^ii, r{y, s) = ^/uj + ^\^n '^^^^ , (6.17a) 

'ikfuj — 1 

cpi = —kicsch^i. (6.17b) 

Similarly, we may apply the method of variation of constant to find the solutions of the 
ACHESCS (5.9) by using the solution (6.17) of (5.3) and corresponding eigenfunction. Taking 
ai in (6.5) to be time-dependent functions ai(s) and requiring that 



Q(y, s) = 2k^csch ^i, r{y, s) = ^/uJ -\ — 2 — , (6.18a) 

4:kfio — 1 

i)! = I3i{s)kicschii (6.18b) 
satisfy the system (j5.9p for = 1, we find that ai(s) can be an arbitrary function of s and 



/3i(.) = Y^^y^M^)- (6.19) 
So the one-cuspon solution of the CHESCS (2.5) with iV = 1 and Ai = A;? - > is obtained 
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with respect to the variables {y,s) from (5.12) 



q{y, s) = ^(1 + 4fc2^_ / ) . (6-20a) 



/ N ^.ifuj^csch^Ci (f. X 

^^2/' - (1 _ 4kfLo){-l + Akju + AkjuJcschHi) ' ^ ' 

^iiy,s)= ^ (6.20c) 
V^(l - 4:kluj){-l + 4A:2w + Akfujcsch'^Ci) 



x{y,s) = ^ + 2 In l-^fci^^coth^ ^Od) 
^yu> 1 + 2Ki-y/ti; coth^i 



In Fig 3, we plot the one-cuspon solution of u, (pi. 




Figure 3. Single cuspon solution for u and the eigenfunction ipi when w — 1, fci = 1, ai{s) — 
-2s, s = 2. 

Similarly, we can apply the method of variation of constant to find the N-cuspon solution of 
(2.5) with Xi = kf - ^ < 0, i = l,--- ,N from (f5l2]l . where 

r{y,s) = V^-2[logWii^u'i>2,--- ,^N)]ysk=^^, (6.21a) 



(1 - 4fe/a;) /n (k - A:,^)^i(^i, *2, • • • , * 



3+^ 



N 



Further more in the same way, we can fnd mixed A:i-soliton-A;2-cuspon solution for (2.5) with 
N = ki + k2, Xi = kf - ^ > 0, i = 1, • • • , fci and Ai = A;2 - ^ < 0, i = Ki + 1,- ■ ■ ,ki + k2, by 
using (6.6) and (5.12). 



6.3 The multipositon solutions 

Let A = —k"^ — Aj = —kf — i = 1, - ■ ■ ,N, and take 



= sin^i, i is an odd number, (6.22a) 
= cos^i, i is an even number. (6.22b) 
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where 



$j = cos^i, i is an odd number, (6.23a) 
$j = sin^i, i is an even number. (6.23b) 



i=i j=i 



— C\k=ki- 



For AT = 1, we have 



*i = sin^i, ^^^^ = 71 cos ^1, (6.24a) 



and 



Wi(^'i, *S'^) = -A;i7i + ^ sin 26, (6.25a) 

1^2(^1, ^r^^i) = -2A;?sin6. (6.25b) 

Then the one-positon solution of (5.3) and the corresponding eigenfunction for (5.1) is given by 
(6.2) with iV = l,mi = 1, 

Q{y, s) = -2[logWi]yy, (6.26a) 

r{y, s) = - 2[logWi]ys. (6.26b) 

V;i(y,s,Ai)=/3i^, (6.26c) 
where ai and /?i are arbitrary constants. 




Figure 4. one-positon solutions for u and the eigenfunction ipi when w = 0.01, ki = 1, ai(s) = 
-2s, s = 2. 
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By using the method of variation of constants, which means we change ai and /3i into ai(s) 
and Pi{s), we obtain the one-positon solution for the CHESCS (2.5) with N = = —k\ — ^ 
from (5.12), where 

f(y, s) = y/u- 2[logWi]ys\^^=^^ , (6.27a) 



ki{l + 4klco)Wi^ 



= 7rTr-^rr2-TT?rl7i=7i' (6.27b) 



li = ai{s)+y+ ^ ^2 ~ ^_L^l,2 • 6.27c) 

(1 + 4/jja;)^ 1 + 4/jja; 

where q:i(s) is an arbitrary function of s. 

In Fig 4, we plot the one-positon solution of u and (pi. 

The positon solution of CHESCS is long-range analogue of soliton and is slowly decreasing, 
oscillating solution In the same way we can find N-positon solution for (2.5). For a detailed 
discussion on positon solution we refer to the reference 



For N, we have 



where 

li 

We find that 



^'j-^^ = 7jCOS^i, i is an odd number, (6.28a) 
^'1^'* = — 7iSin^j, i is an even number. (6.28b) 



Wi = W{^iMi\--- ,*iv,^S^^), (6.29a) 

= W(*i, • • • , *iv, *S^\ $i), (6.29b) 

and the N-positon solution of (5.3) and the corresponding eigenfunction for (5.1) is given by 

Q{y,s) = -2[logWi]yy, (6.30a) 

r{y, s) = - 2[logWi]ys. (6.30b) 
W2 

'ipj{y,s,Xi) = Pi—,i = I,-- - ,N, (6.30c) 

where ai and /3j are arbitrary constants. 

By using the method of variation of constants, we obtain the N-positon solution for the CHESCS 
(2.5) from (5.12), where 

r{y, s) = - '^[logWi]ys\ji=ji, (6.31a) 

^"'(^'^^ = k^ii +4kf.) n,^. (k, + a/ (6.31b) 



7i = 



where ai{s) are arbitrary functions of s. 
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6.4 The multinegaton solutions 

Let X = k'^ - ^ > 0, Xi = kf - ^ > 0, z = 1, • ■ • , AT, and take 

'^i = sinh^i, i is an odd number, (6.32a) 

*i = cosh^i, i is an even number. (6.32b) 

= e«% (6.33) 



where 



■ 1-1 k ki 

1=1 j=i 



then we have 



— C\k=ki J 



*i = sinh^i, '^'i^ = 71 cosh^i, (6.34a) 
Ci = h{y+ f\ ). (6.34b) 

^^ = "^ + ^+(4^:^731)2 (6-34C) 



and 



Wii^i, ^P) = -fci7i + ^ sinh2a, (6.35a) 
W"2(*i,*i^\^i) = 2A;?sinhCi, (6.35b) 

Then the onc-ncgaton solution of (5.3) and the corresponding eigenfunction for (5.1) is given by 

Q{y,s) = -2[logWi]yy, (6.36a) 
r{y, s) = V^- 2[logWi]ys. (6.36b) 

,^i(y,5,Ai)=/3i^, (6.36c) 

where a and /3 are arbitrary constants. 

By using the method of variation of constants, we obtain the one-negaton solution for the 
CHESCS (2.5) with N = l,Xi = kf- ^ from (5.12), where 

f(y, s) = y/uj - 2[logWi]ys\y^=fj^, (6.37a) 



71 = «i(s)+y + 77ixT — T\2+7iX: — 7- {6.37c) 



-16fcfu;5/2g 2u;3/2^ 
(4A;2a; - 1)2 + 4A;2a; - 1 

where a{s) is an arbitrary function of s. 

In Fig 5, we plot the one-negaton solution of u and ipi 
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Figure 5. one-negaton solution for u and the eigenfunction c^i when w = 0.01, fei = 1, a\{s) = 
2s, s = 2. 

Similarly, we can find N-negaton solution for CHESCS (2.5) with Xi = kf-^, i = l,--- ,N. 
For N, we have 



^'j-^^ = 7jCosh,^i, i is an odd number, 
^'•^^ = 7jSinh^j, i is an even number. 



(6.38a) 
(6.38b) 



where 



dk 



+ 



{4kfuj-lf Akfuj-l 



We find that 



(6.39a) 
(6.39b) 

and the N-negaton solution of (5.3) and the corresponding eigenfunction for (5.1) is given by 



Q{y,s) = -2[logWi]yy, 
r{y, s) = ^/u; - 2[logWi]ys. 

W2 

il)j{y,s,\i) = (3i — ,i = ,N, 



(6.40a) 
(6.40b) 

(6.40c) 



where a,; and arc arbitrary constants. 

By using the method of variation of constants, we obtain the N-negaton solution for the 
CHESCS (2.5) from (5.12), where 



r(y, s) = - 2[logWi]ys\^i=^i, 
- 2u: 1 



I7i=7i' 



li = n*^^* ~ kjf^ii^) + y + 



-16fc?a;5/2s 2a;3/2^ 



where 0^(5) are arbitrary functions of s. 

7 Conclusion 



{4kfuj - 1)2 Akfu - 1 



(6.41a) 
(6.41b) 

(6.41c) 



The CHESCS and its Lax representation are derived. Conservation laws are constructed. It is 
reasonable to speculate on the potrential application of CHESCS, that is, CHESCS may describe 
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the interaction of different solitary waves in shallow water. Since SESCS can be regarded as soliton 
equations with non-homogeneous terms, we look for explicit solutions by using the method of 
variation of constants. By considering a reciprocal transformation, which relates CH equation to 
an alternative of ACH equation, we propose a similar reciprocal transformation, which relates the 
CHESCS to ACHESCS. By using the Darboux transformation, one can find the n-soliton and n- 
cuspon solution as well as n-positon and n-negaton solution of alternative ACH equation. Then 
by means of the method of variation of constants, we can obtain N-soliton, N-cuspon, N-positon 
and N-negaton solutions of the ACHESCS. Finally, using the inverse reciprocal transformation, we 
obtain N-soliton, N-cuspon, N-positon and N-negaton solutions of the CHESCS. 
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